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RESEARCH MEMORANDUM 

SUPERSONIC AERODYNAMIC CHARACTERISTICS OF 
A LOW -DRAG AIRCRAFT CONFIGURATION HAVING AN ARROW 
WING OF ASPECT RATIO 1.86 AND A BODY OF 
FINENESS RATIO 20 

By Warren Gillespie, Jr. 

SUMMARY 


A free-flight rocket-propelled model investigation was conducted at 
Mach numbers of 1.2 to 1.9 to determine the longitudinal and lateral aero- 
dynamic characteristics of a low-drag aircraft configuration. The model 
consisted of an aspect-ratio- 1.86 arrow wing with 67 * 5 ° leading-edge sweep 
and NACA 65 AOOL airfoil section, and a triangular vertical tail with 60° 
sweep and NACA 65 AOO 3 section, in combination with a body of fineness 
ratio 20. Aerodynamic data in pitch, yaw, and roll were obtained from 
transient motions induced by small pulse rockets firing at intervals in 
the pitch and yaw directions . 

From the results of this brief aerodynamic investigation, it is 
observed that very slender body shapes can provide increased volumetric 
capacity with little or no increase in zero-lift drag, and that body 
fineness ratios of the order of 20 should be considered in the design of 
long-range supersonic aircraft. The zero-lift drag and the drag-due-to- 
lift parameter of the test configuration varied linearly with Mach number. 
The maximum lift-drag ratio was 7 *0 at a Mach number of 1.25 and decreased 
slightly to a value of 6.6 at a Mach number of 1.8l. The optimum lift 
coefficient, normal-force-curve slope, lateral-force -curve slope, static 
stability in pitch and yaw, time to damp to one-half amplitude in pitch 
and yaw, the sum of the rotary damping derivatives in pitch and also in 
yaw, and the static rolling derivatives all decreased with an increase in 
Mach number. 

Values of certain rolling derivatives were obtained by application 
of the least-squares method to the differential equation of rolling motion. 
A comparison of the experimental and calculated total rolling -moment - 
coefficient variation during transient oscillations of the model indicated 
good agreement when the damping-in-roll contribution was included with the 
static rolling -moment terms . 
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INTRODUCTION 


Several methods have been developed where oy the drag of aircraft 
configurations can be reduced at supersonic speeds. (See refs. 1 to 6 , 
for example.) In general, these methods require either the indentation 
and special contouring of the body in the region of the wing or the 
application of twist and camber to the wing. It may be well to consider 
a simpler approach to the problem of obtaining a low-drag aircraft con- 
figuration suitable for flight at supersonic speeds. For example, con- 
temporary high-speed airplanes have body fineness ratios of the order of 8 . 
The investigation of reference 7 reported in 1951 showed that parabolic 
bodies of fineness ratios 9 to 18 had approximately equal drag at low 
supersonic speeds. However, when based on volume to the two-thirds power 
instead of the usual area reference, the drag coefficient at a Mach number 
of 1.4 for the parabolic bodies was shown to decrease as the fineness ratio 
increased to 25 . 

For the present test a body of fineness ratio 20 was combined with a 
4-percent-thick arrow wing of aspect ratio 1.86. The body was made 
cylindrical in the region of the wing and the overall axial progression 
of total cross-sectional area was moderate. The use of body indentation 
as such was avoided. The purpose of the test was to determine the aero- 
dynamic characteristics of the resulting slender configuration at super- 
sonic speeds and at lifting conditions . The model was flight-tested at 
the Langley Pilotless Aircraft Research Station at Wallops Island, Va. 


SYMBOLS 



normal-force coefficient, 


side -force coefficient, 


W/S 
g H 

ay W/S 
g T 


£Ly W /S 

axial -force coefficient, — — 

g T 

lift coefficient, Cjj cos a + sin a 

drag coefficient, -Cy cos a + Cjj sin a 

minimum drag coefficient (at Cp, = 0) , -Cy at Cp j = 0 
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l/d 

Cm 

c n 

C l 


Cn - 


lift-drag ratio 

pitching -moment coefficient about center of gravity 

yawing-moment coefficient about center of gravity based on 
wing area and span 

rolling -moment coefficient about body center line 


,0a . 


CN=0 


■'“C 


N 


/o c m 

static stability parameter in pitch, 

\ oC 


N/Cn-0 




= 


P VOp 


' Cy=0 


'“c. 


static stability parameter in yaw, 




i) 


\ Cm q. + Cffl d 


K - s*) 


'SCn\ 


sum of rotary damping derivatives in pitch, 

Cm = — and C m . = — — 

m q ^ m a 


2V 


^ac 

2V 


sum of rotary damping derivatives in yaw, 
C„ = — and C„. = 

2V 2V 


C l/2 


time for a transient oscillation to damp to one -half amplitude, 
sec 


period of oscillations, sec 


- 


ac- 


Cl 


p 03 

0c 7 

= 1 

P^Q' 0a 
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c 2 r 


C 


*3 


2V 

Sc 2 





2V 


a n> a y> 


q 

V 

M 

R 

W 

a 

a 

3 

3 

r 

0 

P, 0 
P 

e 

e 


normal, lateral, and axial accelerations, respectively, 
ft/sec 2 

acceleration due to gravity, 32.2 ft/sec 2 
dynamic pressure , lb/ sq ft 
velocity, ft/sec 
Mach number 

Reynolds number based on a length of 1 foot 
weight of model, 111.1 lb 

angle of attack at model center of gravity, deg 
rate of change of angle of attack, radians/sec 
angle of sideslip at model center of gravity, deg 
rate of change of angle of sideslip, radians/sec 
rate of change of flight-path angle , radians/ sec 
angle of roll, deg 
rolling velocity, radians /sec 
rolling acceleration, radians/sec 2 
angle of pitch, deg 

angular velocity in pitch, radians/ sec 


CONFIDENTIAL 


NACA RM L57A25 


CONFIDENTIAL 


5 


8 angular acceleration in pitchy radians/sec^ 

angle of yaw, deg 

r, \jr angular velocity in yaw , radians/sec 

S total wing area to body center line, 4.31 sq ft 

c local wing chord 

c wing mean aerodynamic chord, 2.03 ft 

b wing span, 2.83 ft 

x* chordwise distance back from leading edge of local chord 

y spanwise distance of local chord out from body center line 

Iy moment of inertia in pitch about center of gravity, 

12.2 slug -ft 2 

l£ moment of inertia in yaw about center of gravity, 12.3 slug-ft 2 

Iy moment of inertia in roll about model center line, 

0.314 slug-ft 2 

^•XZ product of inertia, assumed equal to zero 

The positive directions of the angles and coefficients are shown in 
figure 1. 


MODEL 


A drawing of the model is shown in figure 2 and photographs of the 
model are presented in figures 3 and 4. The fuselage ordinates are 
listed in table I, and physical characteristics of the model are listed 
in table II. The configuration for this test consisted essentially of 
an arrow wing of aspect ratio 1.86 with 67.5° leading-edge sweep and 
NACA 65A004 airfoil section attached at body-center-line height to the 
cylindrical midsection of a slender body of fineness ratio 20. The 
model was somewhat similar to the large body configuration, model 5 of 
reference 8. A triangular vertical tail with 60° leading-edge sweep and 
NACA 65AOO3 airfoil section provided directional stability. The tail 
was mounted on top of the body to simulate an airplane conf iguration . 

The ratio of fuselage frontal area to wing plan-form area was O.O32. 
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The angle of incidence of the wing with respect to the body center line 
was zero. The wing dihedral was also zero. 

The model was of metal construction with a solid aluminum-alloy 
wing. Six pulse rockets were carried within the forward and rearward 
fuselage sections, with four firing in the pitch direction and two in 
the yaw direction. The model also carried an eight-channel telemeter 
with angle -of -attack angle-of -sideslip, accelerometer, and rate-of-roll 
instruments. The model was externally boosted by two Deacon rockets. 

An underslung adapter was used to couple the model and booster. A sup- 
port fitting, shown in figure 2, extended below the fuselage and remained 
with the model. 


TEST 


A wing panel and the vertical tail were statically tested to meas- 
ure the streamwise wing twist due to loading concentrated along the 
50-percent-chord line. The flexibility of these model components is 
presented in figures 5 to 7 • 

The model was flight tested at Mach numbers from 1.2 to 1.9 at the 
Langley Pilotless Aircraft Research Station at Wallops Island, Va. Data 
were obtained during ascent of the model after separation from the booster. 
A smoke trail of short duration was generated from a chemical solution con- 
tained in the end of the model which aided in tracking the flight. Aero- 
dynamic data in pitch, yaw, and roll were obtained from transient oscil- 
lating motions induced by pulse rockets firing at intervals in the pitch 
and yaw directions. The telemeter system permitted the measurement of 
angles of attack and sideslip; normal, lateral, and longitudinal accelera- 
tions; angular accelerations in pitch and roll; and rolling velocity. The 
velocity obtained from a CW Doppler radar set (corrected for wind velocity) 
was used in conjunction with tracking radar and radiosonde data to calcu- 
late Mach number, Reynolds number, and dynamic pressure. The variations 
of the free-stream Reynolds number per foot of length and dynamic pressure 
with Mach number are shown in figure 8. Variations of the angle of attack 
with induced sideslip angle caused by pitch pulses are shown in figure 9 . 
Likewise, the variations of the induced angle of attack with sideslip 
angle caused by yaw pulses are shown in figure 10. The variations are 
for the maximum oscillations obtained after a pulse. 


ACCURACY AND CORRECTIONS 


Errors in the absolute value of a telemetered quantity are thought 
to be within ±1 percent of the range of the instrument. At a Mach number 
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of 1.5 the resulting errors in the normal-, lateral-, and axial -force 
coefficients have been calculated to be within ±0.01, ±0.001, and ±0.001, 
respectively. Mach number is estimated to be accurate within ±1 percent 
and dynamic pressure within ±2 percent. Experience in the use of the 
air-flow indicator shows that an error of ±0.3° is probable. 

In order to avoid error, in the determination of the drag polars, 
that could result from either external or internal misalinement of the 
longitudinal (axial) accelerometer instrument when subjected to normal 
acceleration, the angularity of the mounting base in the model was 
measured. The instrument itself was calibrated while subjected to normal 
acceleration. The base of the accelerometer was ground to reduce the 
response of the instrument to normal-force interaction. The residual 
internal instrument error due to normal acceleration and the external 
misalinement of the instrument mounting base were accounted for in the 
data reduction. 

An additional source of inaccuracy in the final results may be the 
induced lateral motions following a pitch pulse or the induced pitch 
motions following a yaw pulse. The relative magnitude of the induced 
lateral motions to pitch motions increased with an increase in Mach number. 
However, cross-coupling effects on the data presented are believed to be 
small . 

Measurements obtained from the flow indicator were corrected for 
pitching and yawing velocities and for flight-path curvature. Position 
corrections were made to measurements obtained from the normal, lateral, 
and longitudinal accelerometers mounted near the center of gravity of 
the model. 


•The instantaneous pitching moment was measured by means of an angular 
accelerometer. The pitching moment due to angle of attack is given by the 
following expression: 


However, for the present test the rotary-damping terms were negligible, 
and the pitching moment due to angle of attack was calculated by the fol- 
lowing simplified expression: 


ANALYSIS 
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The oscillations in pitch resulting from the pitch-pulse rockets . 
have been analyzed assuming two degrees of freedom. A similar analysis 
was made for the oscillations in yaw caused by the yaw -pulse rockets. 
Values of and were calculated using the following expressions: 




57 • 3q.Sc 



These values were divided by corresponding values of CN a and Cyp to 
obtain the static stability parameters and Cn^- Rotary damping 

derivatives were calculated as follows : 



instantaneous rolling moment 'was also measured by means of an 
angular accelerometer. Rolling -moment derivatives were obtained by 
application of the method of least squares to the differential equation 
of rolling motion. Determination of the rolling -moment derivatives is 
explained in the appendix. 


RESULTS AND DISCUSSION 


The aerodynamic test results are presented in figures 11 to 25 for a 
configuration having a wing and a vertical tail with flexibility charac- 
teristics that could be representative of a typical aircraft in this speed 
range. No aeroelastic corrections have been made to the measured data 
obtained during free-f light of the model. 
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Trim 

Figure 11 presents the trim measurements for the model. Because 
the model was not perfectly symmetrical or because of measurement inaccu- 
racies, the trim values for angle of attack, normal-force coefficient, 
angle of sideslip, rolling velocity, and lateral -force coefficient are 
slightly different from zero. The trim angle of attack and normal-force 
coefficient were constant with change in Mach number. The trim angle of 
sideslip, rolling velocity, and lateral -force coefficient all decreased 
with increasing Mach number . 


Drag 

Drag polars were obtained at Mach numbers of 1.25, 1.46, I.69, and 
1.8l and are shown in figure 12. Plots of normal-force coefficient 
against axial-force coefficient are plotted also. The data indicate a 
reduction in axial-force coefficient with increase in normal -force coef- 
ficient. This reduction may be due in part to some suction on the highly 
swept leading edge of the wing, and also to less unfavorable interference 
from the wake of the flow indicator which probably induces a turbulent 
boundary layer well forward on the body of the model, particularly at 
zero angle of attack. In this connection the results of references 9 
and 10 show that the drag at zero lift of a 60° delta-wing — body configu- 
ration (of similar size to the present test model) was 12 to 16 percent 
higher with an air-flow indicator. 


The drag coefficient at zero lift is plotted against Mach number in 
figure 15(a) and is seen to decrease linearly with increase in Mach num- 
ber. A comparison is made with the large body configuration, model 5; of 
reference 8 which had the same wing plan form and maximum wing thickness 
and almost the same ratio of body maximum cross-sectional area to wing 
area. The body fineness ratio was 14.9, however, compared with 20 for 
the present model. The drag at zero lift of the two models is almost the 
same. The present .test model has a considerably larger ratio of 
(Fuselage volume ) 2 ' : 

Wing area 

ratio (hereinafter called the relative fuselage volume) is 0.202 for the 
present test model and 0.148 for the reference model. This drag compari- 
son is somewhat surprising in view of the fact that the present model was 
not aerodynamically "clean" inasmuch as it had six pulse -rocket holes in 
the fuselage in addition to a sting -mounted flow indicator. 


than the model of reference 8. The value of this 


The larger 60° delta-wing — body configuration (model 4 of ref. 9) 
had very nearly the same ratio of body maximum cross-sectional area to 
wing area as the present test model. A direct comparison of the zero- 
lift drag of these two models is made in figure 15(a) and indicates 
slightly lower drag for the larger size model of reference 9- However, 
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if allowance is made for the higher test Reynolds number of model 4 of 
reference 9 and the relative "cleanness" of these two models , then the 
present test model is indicated to have approximately the same drag. 
References 8 and 9 further indicate that the wing-with-interference drag 
of these two models is also approximately equal. The relative fuselage 
volume for the model of this test is, of course, much larger, being 0.202 
for the present model but only 0.113 for model 4 of reference 9. It is 
observed, on the basis of the foregoing comparisons at zero lift, that 
very slender body shapes can provide increased volumetric capacity with 
little or no increase in drag. 

The variation of the drag-due-to-lift parameter dC D /dC L 2 with Mach 

number is linear. (See fig. 13 (b).) This wing plan form is not an opti- 
mum one, particularly at low supersonic speeds. Comparison with the 60° 
delta-wing model of reference 10 shows lower drag-due-to-lift values for 
that model. However, at a Mach number of 1.6 the drag-due-to-lift param- 
eter of the present model is only 4 percent higher and probably would be 
equal at a Mach number of 1.7. 


As a consequence of the linearity of both the variation of the zero- 
lift drag and the variation of the drag-due-to-lift parameter with Mach 
number, the drag of the test model at lift can be represented with good 
accuracy over the test range of Mach number by an expression of the fol- 
lowing form: 


C D = 


( Cd o)m=i . 


dC T 


25 


dM 


^(M - 1.25) 


f dCp 

v^c L 2 


f ac D 

l dC L 2 > 

dM 


(M - 1.25) 


M=1 . 25 


Such an expression might be of value in simplifying the preliminary per- 
formance calculations encountered in the determination of an optimum 
supersonic aircraft with the restriction that the configuration be not too 
far different from that of the present test model for which this result is 
specifically applicable. 


Lift -Drag Ratio 

Figure 14 presents the variation of lift-drag ratio with lift coef- 
ficient obtained at Mach numbers of 1.25, 1-46, I.69, and 1.8l. The 
dashed-line extensions of the plots at the two higher Mach numbers were 

obtained using the expression C D = + ^dC D /dC L 2 ^C L 2 and figure 13. 

The points were plotted using both positive and negative regions of the 
lift-drag data. Maximum lift-drag ratios of 7.0 to 6.6 are indicated to 
occur at an optimum lift coefficient of approximately 0.2. The variations 
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of maximum lift-drag ratio and optimum lift coefficient with Mach number 
are shown in figure 15 . Both quantities decrease nearly linearly with 
increasing Mach number. 

The maximum lift-drag ratio of the present model compares favorably 
with the results obtained for the aspect-ratio- 3.5 swept-wing airplane 
configuration of references 11 and 12. This referenced configuration had 
very nearly the same ratio of maximum body cross-sectional area to wing 
area as the present model, and like the present model had a cylindrical 
fuselage in the region of the wing intersection. The relative fuselage 
volume is less, however, being 0.15 for the fineness-ratio-14.3 fuselage 
of reference 11. It can therefore be stated that body fineness ratios of 
the order of 20 should be considered in the design of long-range supersonic 
aircraft. 


Normal Force and Pitching Moment 

Figures l 6 to 18 present plots of normal-force and pitching -moment 
coefficients and summarize the variations of the normal -force -curve and 
pitching -moment-curve slopes with Mach number. Figure l 6 shows that the 
variation of normal -force coefficient with angle of attack is essentially 
linear for small angles of attack. However, the data for a Mach number 
of 1.25 show that the slope of the curve increases at the higher 

values of Cu and a. The data for Mach numbers of 1.46, 1.69-, and 1.8l 
do not extend far enough in the Ch and a ranges to indicate whether a 
similar increase in Cjj occurs. However, the force data of reference 13 

for a 68.4° delta wing show that an essentially linear variation would be 
expected at a Mach number of 1.9 up to an angle of attack of about 8 ° 
where the slope should begin to decrease because of separation effects. 

The variation of normal -force coefficient with pitching-moment coeffi- 
cient presented in figure 17 is approximately linear over the range of 
the test conditions. The variation of the normal -force -curve slope 

with Mach number shown in figure 18(a) is linear and decreases from a 
value of 0.041 at a Mach number of 1.25 to 0.035 at a Mach number of 1.8l. 
Experimental values of Cjj are approximately 5 percent lower than the 

0G 

values obtained when using the theoretical method of reference 14. This 

comparison indicates very little probable loss in Cm due to wing 

a 

flexibility. A rough estimate based on the aeroelastic analysis of the 

3 -percent-thick, 60 ° delta wing used on the model of reference 10 gives 

a probable reduction of Cm from rigid-wing values of only 4 percent. 

a 

Consequently, a more detailed aeroelastic analysis has not been made for 
the present test model, since the effects of aeroelasticity are probably 
small . 
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The variation of the static stability parameter C^-,^ with Mach 

number is shown in figure 18(b). The experimental results obtained by 
two methods are in close agreement. It is indicated by this agreement 
that lateral oscillations which accompanied the longitudinal motions had 
a negligible effect on the longitudinal period. The result calculated 
using the method of reference 14 compares favorably with the experimental 
curves, but does not show the gradual reduction in static stability as 
Mach number increases. This reduction noted in the tests is probably 
caused by greater wing aeroelasticity for conditions of increased dynamic 
pressure which occurred at the higher test Mach numbers. 


Longitudinal Dynamic Stability 

Figure 19(a) shows that the time for the pitching oscillation to 
damp to one -half amplitude decreased with an increase in Mach number, or 
that the total damping increased with Mach number. One would expect a 
more uniform decrease in j with Mach number rather than the leveling- 

off tendency shown in the figure at the higher Mach numbers. This effect 
is reflected in figure 19(b) which shows negligible rotary damping in 
this region. The theory and experimental tests of references 10, 15, and 
16 indicate that at a Mach number of 1.8 the damping derivatives (0^+0^ 

should have a value of about -0.8 to -0. 5. The slope of the curve of fig- 
ure 19(b) is four times greater than the results of references 10, 15, and 
l6 indicate. The average value of the curve of figure 19(b) is, however, 
in agreement. It should be pointed out that the experimental accuracy of 
the damping derivatives ^C^ + is very poor, because these deriva- 

tives are obtained from the difference of two numbers having the same 
order of magnitude. The important point to be made is that the level of 
the total pitch damping for this tailless (no horizontal tail) configura- 
tion was low, being only one-third that determined for the model of ref- 
erence 12 which had a horizontal tail. 


Side Force and Static Directional Stability 

Plots of side -force coefficient against angle of sideslip are pre- 
sented in figure 20 for Mach numbers of 1.25, 1.46, 1.59, 1.69, 1.8l, and 
1.86. For the small range of the measurements, the variation of Cy 
with p is linear. The slopes obtained from the curves of figure 20 have 
been used to obtain the variation of Cy^ with Mach number shown in fig- 
ure 2l(a). The variation is approximately linear. The static stability 

parameter C n obtained from periods of the yaw pulses is also plotted 

Cy 

against Mach number in figure 2l(b). Comparison with the corresponding 
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data of figure 18(b) indicates that the aerodynamic center in yaw was 
0.2c farther rearward than the aerodynamic center in pitch. 


Directional Dynamic Stability 

Figure 22(a) shows that the time for the yawing oscillations to damp 
to one-half amplitude decreased with an increase in Mach number. The 
rotary yaw damping ^Cn r - decreased slightly with increased Mach num- 

ber. (See fig. 22 (b).) 


Rolling-Moment Derivatives 

Rolling -moment derivatives fCjp at zero angle of attack and C^p 

per degree angle of attack) were obtained by application of a least- 

squares method to the differential equation of rolling motion. The 
method is outlined in the appendix. Although the method is theoretically 
capable of also determining the derivatives C^ and /Cj - C^j, accu " 

rate values of these rotary derivatives could not be determined. Esti- 
mates indicate that the contributions of these terns, particularly of 
^Cj r - , are small in comparison with the contributions of the static 

rolling -moment derivatives ( C 7 | and C 7 to the total rolling 

moment experienced by the model. This is a fortunate situation, and it 
appears that those derivatives which have a greater influence on the 
motion of a particular configuration will be the ones that can be more 
accurately evaluated by this method of data reduction. 


The least-squares method is applicable irrespective of the uniformity 
of the lateral motions. Simultaneously occurring lateral and longitudinal 
(or cross-coupled) motions can be utilized for purposes of stability- 
derivative evaluation. The derivative C^p can be broken down to its 

fundamental parts, (C 7 \ and C 7 . Thus, the motion restrictions 

v Vo-o 

necessary to the proper employment of other methods such as the graphical 
vector method (used in refs. 17 and 18, for example) are greatly relaxed 
or avoided, and the stability derivatives q a 

obtained in lieu of the single derivative Cjp corresponding to some 

average condition of longitudinal trim. 


Figure 23 presents the values of rolling-moment coefficients obtained 
from rolling motions of the model caused either by pitch or yaw pulses . A 
reduction of the absolute values with increase in Mach number is noted. 
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The experimental results are compared with the theoretical variations 
calculated by using the appropriate formulas of references 19 and 20 . . 

The agreement is seen to be generally satisfactory although the theories 
predict somewhat higher values in both cases . For this configuration the 
vertical tail was the largest contributor to (Clp) ^ whereas the wing 

was the largest contributor to Cj R . It should be noted that the theo- 

p ^(X 

retical calculations did not include any interaction effects between com- 
ponents of the configuration tested. Apparently, such effects were small 
for the conditions of the present test. 

Since examination of the transient motions which occurred as a result 
of the yaw pulses showed that amplitude ratios and phase relationships 
could be determined, the vector method of analysis employed in reference l8 
was also used to determine values of C^. The results of this analysis 

are also plotted in figure 23 where a comparison is made with the previ- 
ously determined values of The agreement is good, probably 

because the trim angle of attack was nearly zero. 

The vector analysis also gave values of damping-in-roll parameter 
C, of -0.12 at a Mach number of 1.59 and of -0.14 at a Mach number of 

1.86. These values compare favorably with the level of values obtained 
by the least-squares method and also with the results reported in refer- 
ence 21. 

By using the values of the rolling-moment coefficients obtained from 
the foregoing analysis, comparisons of experimental and calculated rolling- 
moment-coefficient variations with sideslip angle were made. These com- 
parisons are presented in figure 24 at Mach numbers of 1.25, 1.46, 1 . 69 , 
and 1.8l for the case of the model pulsed in pitch; and in figure 25 at 
Mach numbers of 1.59 and 1.86 for the case of the model pulsed in yaw. 

The agreement is generally good when the contributions of 

C-> and C-, are summed. The rolling-moment contribution of the gyro- 

L J3 y CL \ 

scopic reaction ^namely, - ly ®^ J was found to be negligible in the 

determination of the total rolling -moment coefficient Cj. 


SUMMARY OF RESULTS 


Results obtained from a flight test of a low-drag aircraft configu- 
ration at supersonic speeds lead to the following observations : 
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1. Very slender body shapes can provide increased volumetric capacity 
with little or no increase in zero-lift drag. Body fineness ratios of the 
order of 20 should be considered in the design of long-range supersonic 
aircraft . 

2. Maximum lift-drag ratios of 7*0 and 6.6 at Mach numbers of 1.25 
and 1.8l, respectively, were obtained. 

5. The optimum lift coefficient, normal-force-curve slope, lateral- 
force -curve slope, static stability in pitch and in yaw, time to damp to 
one -half amplitude in pitch and in yaw, the sum of the rotary damping 
terms, and the Static rolling derivatives all decreased with, an increase 
in Mach number. 

4. Comparison of the experimental and calculated variation of the 
total rolling-moment coefficient during transient oscillations of the 
model indicated good agreement when the damping-in-roll contribution was 
included with the static rolling -moment terms. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 7, 1957* 
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APPENDIX 

DETERMINATION OF ROLLING DERIVATIVES BY THE 
LEAST-SQUARES METHOD 


In order to utilize the transient rolling measurements obtained 
immediately following the pitch disturbances for the purpose of deter- 
mining rolling derivatives, the least-squares method of data reduction 
was applied to the differential equation of rolling motion. The least- 
squares method is outlined in reference 22, pages 371 and 372. Data from 
both pitch and yaw pulses were analyzed to obtain values of the rolling 
derivatives . 


The total net aerodynamic rolling-moment coefficient at any instant 
during free oscillation is given as follows: 


C Z 


I X P + (ig ~ lyj^V - Ixz^ P® ) 
qSb 


( 1 ) 


For the present model the product of inertia was assumed to be equal to 
zero, and the contribution of the gyroscopic reaction term was found to 
be negligible . The net aerodynamic rolling -moment coefficient was then 
obtained from the following simplified expression: 



( 2 ) 


This net aerodynamic coefficient was next assumed to result from a simple 
addition of particular rolling -moment coefficients . Thus , 


C z = + CKj\|r + CK j| p 


(3) 
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where 


Ki = 



K 3 - (% - C Z .) 

k 4 = % 


c 


jb_ 

2V 


From the telemeter, radar, and radiosonde measurements, sets of data 
were obtained consisting of p, q, (3, a, \jr, p, and V at selected 
times over approximately 1 oscillation corresponding to an average Mach 
number. The largest amplitude oscillations immediately following a pitch 
or yaw pulse were used. Trim conditions for the telemetered quantities 
were determined, and the sets of data corrected to incremental variations 
from trim. The corresponding values of were calculated by using 

equation (2). The following equations can then be written: 


c z = + 

C z = + K2OC232 + ^ 2 ^‘ 3^ f 2 + 


00 


c Z n = K l3n + K&rfin + C n K 3 t n + C n K 4 p n 

/ 

The unknowns are the K f s and the subscripts (the K ! s excepted) 
refer to particular sets of data. Choose as the best approximation to 
the unknowns those values which minimize the sum of the squares of the 
deviations of the observed values from the corresponding values which 
the observed quantity would have if computed from the chosen values of 
the unknowns . The following expression can then be minimized by equating 
to zero the four partial derivatives with respect to Kq_, Kj, 

and Ki|.: 


CONFIDENTIAL 


18 


CONFIDENTIAL 


NACA EM L57A25 



[ c lj - ( K lPj + * C J K 5*J + C J K ‘t p j) 


(5) 


There results the following set of equations 


where 


Ef = + BjIG* 

E^ = 

Ej = AjK-^ + B^K^ 

E||_ = A^K^ + B^K^ 


+ Ci Kj + D-j_K|^ 
+ 

+ C^Kj + D^Kj^ 
+ C^Kj + D4K4 



( 6 ) 



cp P 


D^ 



Ca(3p 







d a 


y (cp) 2 



Equations (6) must be solved simultaneously for the K's and the 
corresponding aerodynamic parameters ( c Zp)o^O> c Zp, a > c 2p> and 

~ C ^p)‘ The accuracy of determination of these parameters will 

depend on the accuracy and extent of the basic measurements and the 
relative importance of the various terms to the rolling motion of the 
configuration under consideration. In the present case values for 
(Ci \ and for C 7 were determined, but only the order of magnitude 

{ bja^O f \ 

of Cj could be determined. The contribution of (Cz r - C^j was esti- 
mated to be negligible, and accurate values for the sum of these two 
damping derivatives could not be determined. 
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TABLE I.- FUSELAGE ORDINATES 


Station, in. 

Body radius, in. 

0 

0 

.67 

.22 

1-33 

.38 

1.67 

.44 

2.33 

.57 

3.33 

.73 

5.00 

.98 

6.67 

1.19 

10.00 

1.54 

13.33 

1.82 

16.67 

2.06 

20.00 

2.23 

22.75 

2.35 

23.33 

2.37 

26.67 

2.45 

30.00 

2.50 

Constant radius 

Constant radius 

63.38 

2.50 

67.43 

2.45 

71.49 

2.37 

75.5^ 

2.23 

79.60 

2.06 

83.65 

1.82 

87.71 

1.54 

91.76 

1.19 

93.79 

.98 

95.82 

•73 

97.04 

• 57 

97.85 

.44 

98.25 

.38 

99.06 

.21 

99.87 

0 
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TABLE II.- CHARACTERISTICS OF MODEL 

Wing: 

Span, ft 2.83 

Area, sq ft L .31 

Aspect ratio 1.86 

Taper ratio 0 

Sweepback of leading edge, deg 67.5 

Sweepback of trailing edge, deg 15 

Mean aerodynamic chord, c, ft 2.03 

Airfoil section NACA 65 AOO 4 

Incidence, deg 0 

Dihedral, deg 

Body: 

Maximum diameter, ft 0.42 

Length, ft 8.32 

Fineness ratio 20 

Vertical tail: 

Span, ft 0.97 

Taper ratio 

Sweepback of leading edge, deg 60 

Sweepback of trailing edge, deg 15 

Airfoil section NACA 65 AOO 3 

Model weight, lb 111.1 

O 

Moments of inertia, slug-ft^ 

In pitch 12.2 

In yaw 12.3 

In roll 0.314 

Center of gravity, percent c behind leading edge of 

mean aerodynamic chord 0.197 
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Wind 




Figure 1.- System of axes. Arrows indicate positive directions; origin 

is at center of gravity. 
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Figure 2 .- Test configuration. All linear dimensions are in inches 
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Figure 3-- Model used in investigation. 
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Figure 4.- Model and booster. 
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Figure 5 .- Wing static deflection resulting from a concentrated load 
applied along 50 -percent- chord line. 
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Load applied at station 




Figure 5.- Concluded. 
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(a) Average twist from 0 to 50 percent local chord. 

Figure 6.- Wing streamwise twist resulting from a concentrated load 
applied along 50-percent-chord line. 


CONFIDENTIAL 


50 


CONFIDENTIAL 


NACA RM L5TA25 


6xlO- 2 



(b) Average twist from 50 to 80 percent local chord. 


Figure 6.- Continued. 
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(c) Average twist from 80 to 100 percent local chord. 
Figure 6.- Concluded. 
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Figure 7.- Tail streamwise average twist resulting from a concentrated 
load applied along 50-percent- chord line. 
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Figure 8.- Flight test conditions. 
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£, deg 

(a) Mach number, 1.24. 

Figure 9 .- Variation of angle of attack with sideslip angle. Model 

pulsed in pitch. 
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(b) Mach number, 1.46. 
Figure 9-- Continued. 
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(c) Mach number, 1.69* 
Figure 9.- Continued. 
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/ 3 , deg 

(d) Mach number, 1.8l. 
Figure 9.- Concluded. 


CONFIDENTIAL 


CONFIDENTIAL 


2 


V>4 

00 



-2 1 1 ^ ^ ^ ^ ^ ^ 1 — i — ■ — ■ - 

-3-2-1 0 I 2 3 

/3, deg 


t 



(a) Mach number, 1.59- 

Figure 10.- Variation of angle of attack with sideslip angle. Model 

pulsed in yaw. 
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(b) Mach number, 1.86. 
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(a) Longitudinal. 
Figure 11.- Model trim. 
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(b) Lateral. 
Figure 11.- Concluded. 


CONFIDENTIAL 


k2 


CONFIDENTIAL 


MCA RM L57A25 



M= 1.25 1.46 1.69 1.81 


- or 


Figure 12.- Variation of normal-force coefficient with axial-force -coeffi- 
cient and lift coefficient with drag coefficient at constant Mach numbers. 


CONFIDENTIAL 


NACA RM L57A25 


CONFIDENTIAL 


43 


c D °t 
C, =0 
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.01 
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Present test 

Ref. 8, model 5 

Ref. 9 , model 4 
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(a) Drag coefficient at C T = 0. 




M 


(b) Drag-due-to-lift parameter dCj^jdCj^ . 

Figure 13.- Drag coefficient at zero lift and drag-due-to-lift parameter 

plotted against Mach number. 
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Figure 14.- Lift-drag ratio plotted against lift coefficient. 
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(a) Maximum lift-drag ratio. 


.3 

.2 

C L, OPT 


0 

1.2 1.4 1.6 1.8 2.0 
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(b) Lift coefficient at (L/D) 

max 

Figure 15.- Maximum lift-drag ratio and optimum lift coefficient plotted 

against Mach number. 
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Figure l6.- Variation of normal-force coefficient with angle of attack. 
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Figure 17 .- Variation of normal-force coefficient with pitching-moment 

coefficient . 
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(a) Normal-force -curve slope 
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(b) Static stability parameter in pitch 

C N 

Figure l8.- Variation of lift effectiveness and static stability parameters 

with Mach number. 
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(a) Time for pitching oscillation to damp to one-half amplitude. 
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(b) Rotary damping derivatives. 
Figure 19 *- Damping in pitch. 
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Figure 20.- Variation of side-force coefficient with angle of sideslip. 


CONFIDENTIAL NACA RM L57A25 


NACA RM L57A25 


CONFIDENTIAL 


51 


-.010 


-008 


C Y -006 

/ 9 

-.004 


-.002 


0 

1-2 1.4 1.6 1.8 2.0 

M 



C 


Y 


(a) Side-force-curve slope Cyp« 
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(b) Static stability parameter in yaw C n 

C Y 

Figure 21.- Variation of side-force effectiveness and static directional 
stability parameters with Mach number. 
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(a) Time for yawing oscillation to damp to one-half amplitude. 
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(b) Rotary damping derivatives. 
Figure 22.- Damping in yaw. 
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Figure 2}.- Values of rolling-moment coefficients obtained by applica- 
tion of least-squares method to the rolling-moment equation. 
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(a) Mach number, 1.25. 

Figure 2k. - Comparison of experimental and calculated rolling-moment- 
coefficient variation with sideslip angle. Model pulsed in pitch. 
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(b) Mach number, 1.46. 
Figure 24.- Continued. 
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(c) Mach number, I.69. 

Figure 2b.- Continued. 
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(d) Mach number, 1.8l. 
Figure 24.- Concluded. 
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(a) Mach number, 1.59» 

Figure 25 .- Comparison of experimental and calculated rolling-moment- 
coefficient variation with sideslip angle. Model pulsed in yaw. 
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(b) Mach number, 1.86. 
Figure 25.- Concluded. 
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